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The axino and the gravitino are extremely weakly interacting candidates for the lightest supersymmetric particle 
(LSP). We demonstrate that either of them could provide the right amount of cold dark matter. Assuming that 
a charged slepton is the next-to-lightest supersymmetric particle (NLSP), we discuss how NLSP decays into the 
axino/gravitino LSP can provide evidence for axino/gravitino dark matter at future colliders. We show that these 
NLSP decays will allow us to estimate the value of the Peccei— Quinn scale and the axino mass if the axino is the LSP. 
In the case of the gravitino LSP, we illustrate that the gravitino mass can be determined. This is crucial for insights 



' into the mechanism of supersymmetry breaking and can lead to a microscopic measurement of the Planck scale. 

1. INTRODUCTION 
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I A key problem in cosmology is the understanding of the nature of cold dark matter. In supersymmetric extensions 

• of the Standard Model, the lightest supersymmetric particle (LSP) is stable if i?-parity is conserved 0]. An electrically 

. and color neutral LSP thus appears as a compelling solution to the dark matter problem. The lightest neutralino 

[ is such an LSP candidate from the minimal supersymmetric standard model (MSSM). Here we consider two well- 

' motivated alternative LSP candidates beyond the MSSM: the axino and the gravitino. 

In the following we introduce the axino and the gravitino. We review that axinos/gravitinos from thermal pro- 

■""...^ ' duction in the early Universe can provide the right amount of cold dark matter depending on the value of the 

'" ^ , reheating temperature after inflation and the axino/gravitino mass. For scenarios in which a charged slepton is the 

I ' next-to-lightest supersymmetric particle (NLSP), we discuss signatures of axinos and gravitinos at future colliders. 

D . 

; 

:^ 2. AXINOS AND GRAVITINOS 

X. 



The axino a 0,0) 1| appears (as the spin-1/2 superpartner of the axion) once the MSSM is extended with the Peccei- 
Quinn mechanism a| in order to solve the strong CP problem. Depending on the model and the supersymmetry 
(SUSY) breaking scheme, the axino mass m~ can range between the eV and the GeV scale The axino is a 

singlet with respect to the gauge groups of the Standard Model. It interacts extremely weakly as its interactions are 
suppressed by the Peccei-Quinn scale 0, /a ^ 5 x 10^ GeV. The detailed form of the axino interactions depends 
on the axion model under consideration. We focus on hadronic, or KSVZ, axion models in a SUSY setting, in 
which the axino couples to the MSSM particles only indirectly through loops of additional heavy KSVZ (s)quarks. 

The gravitino G appears (as the spin-3/2 superpartner of the graviton) once SUSY is promoted from a global to 
a local symmetry leading to supergravity (SUGRA) The gravitino mass m~ depends stro ngly on the SUSY- 

breaking scheme and can range from the eV scale to scales beyond the TeV region 

111 lia IM In particular, 

in gauge- mediated SUSY breaking schemes 0|, the mass of the gravitino is typically less than 100 MeV, while 
in gravity-mediated schemes Q it is expected to be in the GeV to TeV range. Also, the gravitino is a singlet 
with respect to the gauge groups of the Standard Model. Its interactions — given by the SUGRA Lagrangian — are 
suppressed by the (reduced) Planck scale Mpi = 2.4 x 10^^ GeV. Once SUSY is broken, the extremely weak 
gravitino interactions are enhanced (for small values of the gravitino mass) through the super-Higgs mechanism. 
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Figure 1: Relic densities of axino LSPs (left) and gravitino LSPs (right) from thermal production in the early Universe. 

3. AXINOS AND GRAVITINOS AS COLD DARK MATTER 

Because of their extremely weak interactions, the temperature Td at which axinos/gravitinos decouple from the 
thermal plasma in the early Universe is very high. For example, an axino decoupling temperature of Td « 10^ GeV is 
obtained for fa = 10^^ GeV. Gravitinos with rriQ > lOMeV decouple at similar or even higher temperatures. Below 
the decoupling temperature, axinos/gravitinos can be produced in thermal reactions in the hot MSSM plasma. 
The thermal axino/gravitino production rate at high temperatures can be computed in a gauge- invariant way with 
the Braaten-Yuan prescription and hard thermal loop resummation, which takes into account Debye screening in 
the plasma. Assuming that inflation has diluted away any primordial axino/gravitino abundance, axino/gravitino- 
disappearance processes are negligible for a reheating temperature Tr sufficiently below Td. The corresponding 
Boltzmann equation can be solved analytically. This leads to the following results for the relic densities fl^^^h^ of 
stable LSP axinos/gravitinos from thermal production [ifil 



n,h- = 5.5 5^h. (1) 




n^h^ = 0.12 5^1np^ 1 + -^ — ^ (2) 



where h parametrizes the Hubble constant Hq = 100 /i km/s/Mpc, g is the strong coupling, N the number of heavy 
KSVZ (s)quark flavors, and rag the gluino mass. Note that fl^h^ increases with decreasing m— due to enhanced 
gravitino interactions. The axino couplings do not depend on so that flah^ increases with increasing ma- 
in Fig. n we illustrate fig^^/i^ as a function of Tji for different values of rn^^^, where fa/N — 10^"'^ GeV and 
fng{Mz) = 700 GeV. The running of the strong coupling and the gluino mass is taken into account by replacing g and 



m 



g in (CEl) respectively with g{TR) = [g~^{Mz) + 2,HTR/Mz)/{S^^)]-^'^ and mg(Tfl) = [g{TR) / g{Mz)Ymg{Mz). 



where Mz is the Z-boson mass and g"^ {M z) / {'Iit) 0.118. For Tr above Td w 10^ GeV, fig/i^ is given by the 
equilibrium number density of a rclativistic Majorana fermion and thus independent of Tr as shown for rria — 1 keV. 
There will be a smooth transition instead of a kink once the axino-disappearance processes are taken into account. The 
grey band indicates the WMAP result Q on the cold dark matter density {2a error) Q^^^K^ = O.llSlg p^g. Axinos 
give the right amount of cold dark matter for {ma,TR) — (100 keV, 10^ GeV). Higher values of Tr are problematic 
as rua becomes too light to explain structure formation. In contrast, gravitinos could provide the right amount of 
cold dark matter for combinations with a very high reheating temperature such as (m^, Tr) = (100 GeV, 10^" GeV). 
Nevertheless, other gravitino cold dark matter scenarios — such as {ttiq^Tr) = (lOMeV, 10^ GeV) — are also viable. 
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Figure 2: The 2-body decay tr r + a. 



4. COLLIDER SIGNATURES OF AXINOS AND GRAVITINOS 

As a consequence of the extremely weak couplings, the direct production of axino/gravitino LSPs at colliders is 
strongly suppressed. Furthermore^the NLSP typically has a long lifetime, which is subject to cosmological constraints 
as discussed, for example, in []j|,|2(j|. At future colliders one expects a large sample of such quasi-stable NLSPs if the 
NLSP belongs to the MSSM spectrum. Each NLSP will eventually decay into the axino/gravitino LSP. These decays 
can provide signatures of axinos/gravitinos and other insights into physics beyond the MSSM (cf. [21I E3. E3 . and 
references therein). We concentrate on results extracted from Ref. i24ij where more details can be found. 

A significant fraction of the NLSP decays will take place outside the detector and will thus escape detection. For 
the charged slepton NLSP scenario, however, two recent works have proposed how such NLSPs could be stopped and 
collected for an analysis of their decays. It was found that up to O(IO'^-IO^) and C'(IO'^-IO^) of charged NLSPs can 
be trapped per year at the Large Hadron Collider (LHC) and the International Linear Collider (ILC), respectively, 
by placing 1-10 kt of massive additional material around planned collider detectors 

To be specific, we focus on the case where the pure 'right-handed' stau tr is the NLSP. We assume for simplicity 
that the ncutralino-stau coupling is dominated by the bino coupling and that the lightest neutralino is a pure bino. 



4.1. Probing the Peccei-Quinn Scale and the Axino Mass 

In the axino LSP scenario, the total decay rate of the stau NLSP is dominated by the 2-body decay r — > r + a. In 
Fig. [3 we show the corresponding Feynman diagrams for the considered hadronic (KSVZ) axion models. The heavy 
KSVZ (s)quark loops are indicated as effective vertices by thick dots. The decay rate was estimated as |3 

^r.^.e (: - g) (^) (12^)' . ,3, 

where is the mass of the bino and rrif is the mass of the stau NLSP, i.e. ma < rrif < nig. The KSVZ-model 
dependence is expressed by CaVY — 0{1) and the uncertainty of the estimate is absorbed into ^ ~ Thus, from 

the lifetime of the stau NLSP, Tf « 1/T(ty{. to.), an estimate of the Peccei-Quinn scale fa can be obtained [3| 

f. ^ (^) ec;^. (1 - ^ (^) (lO"GeV)' . ,4, 

Indeed, we expect that ruf and rrig will already be known from other processes when the stau NLSP decays are 
analyzed. The dependence on rria is negligible for rria/rnf ^ 0.1. For larger values of m^, the stau NLSP decays can 
be used to determine the axino mass kinematically, i.e., from a measurement of the energy of the emitted tau E-r, 



m? + — 2mfEr , (5) 

where the error is given by the experimental uncertainty on rrif and E-r ■ The determination of both the Peccei-Quinn 
scale fa and the axino mass is crucial for insights into the cosmological relevance of the axino LSP. 
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Figure 3: The 3-body decays tr ^ r + 7 + a (left) and tr ^ r + 7 + G (right). 

4.2. Measuring the Gravitino Mass and the Planck Scale 

In the gravitino LSP scenario, the main decay mode of the stau NLSP is the 2-body decay r r + G. Neglecting 
the T mass, the following tree-level result for the decay rate is obtained from the SUGRA Lagrangian 

with the value of the reduced Planck mass Mpi = (SttGn)"^/^ = 2.435 x 10^* GeV as given by macroscopic mea- 
surements of Newton's constant ^3 G'n — 6.709 x 10"'^^ GeV~^. Thus, the gravitino mass m~ can be determined 
once the stau NLSP lifetime, Tf « l/r(rR —trG), and rrif are measured. This will be crucial for insights into 
the SUSY breaking mechanism. If the gravitino mass can be determined independently from the kinematics via 
rriQ — {mf + m'^ — 2TOf i?,-)^/^, the lifetime Tf can also be used for a microscopic measurement of the Planck scale ^2^ 

Ml,^^^(l-^X . (7) 
If consistent with macroscopic measurements, this would provide evidence for the existence of SUGRA in nature. 

4.3. Distinguishing between Axinos and Gravitinos 

A question arises as to whether one can distinguish between the axino LSP and the gravitino LSP scenarios. For 
nif — 100 GeV and rrig — 110 GeV, for example, the stau NLSP lifetime in the axino LSP scenario can range from 
0(0.01 s) for /a = 5 X 10^ GeV to 0(10 h) for /a = 5 x lO^^ Qgy ^^e gravitino LSP case, the corresponding 
lifetime can vary over an even wider range, e.g., from 6 x 10~^ s for m~ = 1 keV to 15 years for m~ = 50 GeV. Thus, 
both a very short lifetime, Tf ^ ms, and a very long one, Tf ^ days, will point to the gravitino LSP. On the other 
hand, if the LSP mass cannot be measured and the lifetime of the stau NLSP is within the range 0(0.01 s)-O(10 h), 
it will be very difficult to distinguish between the axino LSP and the gravitino LSP from the lifetime Tf alone. 

The situation is considerably improved when one considers the 3-body decays tr — > t -f 7 -|- a/G. The Feynman 
diagrams of the dominant contributions are shown in Fig. O where thick dots represent heavy KSVZ (s)quark loops 
and shaded triangles the set of diagrams given in Fig. |21 From the corresponding differential rates [2J| , one obtains 
the differential distributions of the visible decay products. These are illustrated in Fig. ^ in terms of the quantity 

1 rf^r(?R^T7») 

^cut ^cut\ j_ j„„„a ' « — [S) 



r(TR ^ T7i ; a;^"', Xg"') dxjdcos9 
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Axino LSP Scenario 



Gravitino LSP Scenario 




Figure 4: DifFerential distributions of the visible decay products in t-r. —> t + "f + a (left) and tr ^ r + 7 + G (right). 



where x-y = 2Ej/mf is the scaled photon energy, 6 is the opening angle between the photon and tau directions, and 

r(?R^r7i;x™*,x^"*)= / dx-, dcose— f^- with A, = ^ (9) 

j^jcut J —1 dx^ucosu 

is the respective integrated 3-body decay rate with the cuts > a;™* and cos 9 < 1 — x™*. Note that the quantity (jSJ 
is independent of the 2-body decay, the total NLSP decay rate, and the Peccei-Quinn/Planck scale. 

The figure shows the normalized differential distributions ^ for the axino LSP with ml/m^ < 1 (left) and the 
gravitino LSP with m~ = lOMeV (right), where nif = 100 GeV, = 110 GeV, and x™* = xg"* = 0.1. The contour 
lines represent the values 0.2, 0.4, 0.6, 0.8, and 1.0, where the darker shading implies a higher number of events. In 
the case of the gravitino LSP, the events are peaked only in the region where the photons are soft and the photon 
and the tau are emitted with a small opening angle (0 ~ 0). In contrast, in the axino LSP scenario, the events are 
also peaked in the region where the photon energy is large and the photon and the tau are emitted back-to-back 
(0 ~ tt). Thus, if the observed number of events peaks in both regions, there is strong evidence for the axino LSP 
and against the gravitino LSP. 

To be specific, with 10^ analyzed stau NLSP dec ays, we expect about 165±13 (stat.) events for the axino LSP and 
about lOOzblO (stat.) events for the gravitino LSP |2J|, which will be distributed over the corresponding (x^,cos0)- 
planes shown in Fig.^] In particular, in the region of x-y > 0.8 and cos 6 < —0.3, we expect about 28% of the 165±13 
(stat.) events in the axino LSP case and about 1% of the lOOilO (stat.) events in the gravitino LSP case. These 
numbers illustrate that 0(10'*) of analyzed stau NLSP decays could be sufhcient for the distinction based on the 
differential distributions. To establish the feasibility of this distinction, a dedicated study taking into account the 
details of the detectors and the additional massive material will be crucial, which we leave for future studies. 

Some comments are in order. The differences between the two scenarios shown in Fig. ^become smaller for larger 
values of rrij^/mf. This ratio, however, remains close to unity for the stau NLSP in unified models. Furthermore, 
if rriQ < rria < rrif and T(t —>■ a X) ^ T(t —>■ GX), one would still find the distribution shown in the left panel of 
Fig. ^ The axino would then eventually decay into the gravitino LSP and the axion. Conversely, the distribution 
shown in the right panel of Fig. ^ would be obtained if rria < < rrif and T(t a X) <C T{t — > G X). Then it 
would be the gravitino that would eventually decay into the axino LSP and the axion. Barring these caveats, the 
signatures shown in Fig. ^ will provide a clear distinction between the axino LSP and the gravitino LSP scenarios. 
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5. CONCLUSIONS 



Axino/gravitino LSPs from thermal production in the early Universe can provide the right amount of cold dark 
matter depending on their mass and the reheating temperature after inflation. If a charged slepton is the NLSP, 
future colliders can provide signatures of axino/gravitino LSPs and other insights into physics beyond the MSSM. 
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